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Abstract 

Background: Since the first outbreak in Indonesia in 1926, Newcastle disease has become one of the most 
common and contagious bird diseases throughout the world. To date, enhancing host antibody response by 
vaccination remains the most efficient strategy to control outbreaks of Newcastle disease. Antibody response plays 
an important role in host resistance to Newcastle disease, and selection for antibody response can effectively 
improve disease resistance in chickens. However, the molecular basis of the variation in antibody response to 
Newcastle disease virus (NDV) is not clear. The aim of this study was to detect genes modulating antibody response 
to NDV by a genome-wide association study (GWAS) in chickens. 

Results: To identify genes or chromosomal regions associated with antibody response to NDV after immunization, 
a GWAS was performed using 39,833 SNP markers in a chicken F 2 resource population derived from a cross 
between two broiler lines that differed in their resistance. Two SNP effects reached 5% Bonferroni genome-wide 
significance (P<1.26x10~ 6 ). These two SNPs, rs15354805 and rs15355555, were both on chicken {Gallus gallus) 
chromosome 1 and spanned approximately 600 Kb, from 100.4 Mb to 101.0 Mb. Rsl 5354805 is in intron 7 of the 
chicken Roundabout, axon guidance receptor, homolog 2 {ROB02) gene, and rsl 5355555 is located about 243 Kb 
upstream of ROB02. Rsl 5354805 explained 5% of the phenotypic variation in antibody response to NDV, post 
immunization, in chickens. Rsl 5355555 had a similar effect as rsl 5354805 because of its linkage disequilibrium with 
rsl 5354805 (r 2 =0.98). 

Conclusion: The region at about 100 Mb from the proximal end of chicken chromosome 1, including the ROBOl 
and ROB02 genes, has a strong effect on the antibody response to the NDV in chickens. This study paves the way 
for further research on the host immune response to NDV. 
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Background 

Newcastle disease is a highly contagious viral disease of 
birds affecting over 250 domestic and wild avian species. 
Since Newcastle disease was first discovered in Indonesia 
in 1926, it has spread throughout the world [1]. Newcas- 
tle disease is one of the most serious problems affecting 
poultry industries in many countries because its conta- 
giousness and devastating effects [2-5]. Newcastle disease 
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virus (NDV) has been classified into velogenic, mesogenic 
and lentogenic strains on the basis of their pathogenesis 
and virulence. Chickens with Newcastle disease have severe 
neurological and respiratory signs and show decreased egg 
quality and production [1,6,7]. There is no effective treat- 
ment for Newcastle disease; however, the use of prophylac- 
tic vaccines and maintenance of strict biosecurity measures 
can reduce the likelihood of outbreaks. Thus, the ability of 
chickens to mount an antibody response to NDV plays a 
key role in controlling Newcastle disease outbreaks, why 
understanding the molecular basis of immune response to 
NDV is important for the control of avian Newcastle 
disease. The antibody response to the same virus differs 
between chicken breeds [8], and selection for an antibody 
response may improve disease resistance in chickens [9]. 
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However, there have been few studies on controlling 
Newcastle disease outbreaks from the perspective of 
host resistance. Ten quantitative trait loci (QTLs) have 
been reported to be linked to the antibody response to 
NDV on chicken (Gallus gallus) chromosomes (GGA) 
2, 3, 4, 5, 9, 13, 16, 18, 22 and Z [10,11], but only a few 
causative genes have been identified because of low map 
resolution. The chicken 60K single nucleotide poly- 
morphism (SNP) chip provides an opportunity to in- 
crease precision of QTL mapping by increasing genome 
coverage and map resolution [12]. Using this chicken 
60K SNP chip, some QTLs and candidate genes affect- 
ing chicken growth have been identified at 172-175 Mb 
on GG Al [13] and 71.6-80.2 Mb on GGA4 [14]. 

The aim of this study was to identify QTLs for the 
antibody response to NDV in chickens. A genome-wide 
association study (GWAS) was performed using the 
chicken 60K SNP chip, based on the data from an F 2 re- 
source population that was derived from the cross be- 
tween two broiler lines with differing resistances. 

Results 

Characterization of the antibody response to NDV in the 
F 2 population 

The antibody response to NDV was positive in all 511 
F 2 birds at 41 days after the second immunization, as 
assessed by an enzyme linked immunosorbent assay 
(ELISA) test, and were highly variable. The mean ± 
standard error of the adjusted antibody responses to 
NDV was 3.38 ± 0.05, its coefficient of variance was 
34.0%, and its distribution fitted a normal distribution 
(P=0.5527, Additional file 1: Figure SI). The heritability 
of the secondary antibody response to NDV was esti- 
mated as 0.24 ± 0.08. 

Genome-wide association analysis 

According to the quality control criteria, 505 out of 511 
birds in the F 2 population and 39833 out of 57636 SNPs 
on the Illumina 60K Chicken SNP Beadchip were eligible 
for inclusion in the genome-wide association analysis. 
The informative SNPs in this GWAS were distributed 
across the chicken genome, at an average interval of 
26.4 Kb (Additional file 2: Table SI). The global view of 
P values for all SNPs affecting the antibody response to 
NDV showed that a region on GGA1 was strongly asso- 
ciated with the chicken antibody response to NDV at 41 
days after the second immunization (Figure 1). Two 
SNP effects reached 5% Bonferroni genome-wide signifi- 
cance (P<1.26xl0~ 6 ): rsl5354805 at 100399530 bp of 
GGA1 and rslS3SSSSS at 100994585 bp of GGA1. SNPs 
rsl5354805 and rsl5355555 explained 5% and 4% of the 
phenotypic variation in antibody response to NDV at 41 
days after the second immunization in the F 2 population, 
respectively (Table 1). Fourteen SNPs with suggestive 



association (P<2.51xl0~ 5 ) covered a region of approxi- 
mately 12.2 Mb from 92.2 Mb to 104.4 Mb on GGA1. All 
the SNPs with P values smaller than l.OOxlO" 4 were con- 
centrated in the 89.7-111.0 Mb region of GG Al, except 
for one SNP on GGA10 (Additional file 3: Table S2). 
Collectively, all the SNPs with P<1.00xl0~ 4 explained 8% 
of the phenotypic variation in the antibody response to 
NDV, whereas just two SNPs, rsl5354805 on GGA1 and 
rsl4008095 on GGA10, accounted for 7% of this pheno- 
typic variation. 

Evaluation of the candidate region comprising rsl5354805 
and rsl 5355555 

Strong linkage disequilibrium (LD) in GGA1 occurred be- 
tween SNP markers within a distance of approximately 
200 Kb, and the LD declined to the background level 
(r 2 =0.1) at a distance of about 2 Mb in the F 2 population 
(Additional file 4: Figure S2). As shown in Figure 2, there 
are 12 haplotype blocks in the 99-101 Mb region of 
GGA1, including rsl 5354805 and rsl 5355 555. This region 
contains three protein coding genes: LOC769250, Round- 
about, axon guidance receptor, homolog 1 (ROBOl) and 
Roundabout, axon guidance receptor, homolog 2 (ROB02). 
Rsl5354805 and three additional SNPs (GGAluGA033635, 
rsl3908604 and rsl5354828) formed an LD block at a dis- 
tance of 92 Kb in the F 2 population. All four SNPs are in 
ROB02, and rsl5354805 is located in intron 7 of chicken 
ROB02. Although rsl5354805 and rsl5355555 are in an 
approximately 600 Kb interval, and are not in the same 
haplotype block, the pairwise LD, as represented by the r- 
square value between rsl5354805 and rsl5355555, reached 
0.98. SNP rsl5355555 is located about 243 Kb upstream of 
ROB02. This result suggested that the same underlying 
genetic factor caused rsl5354805 and rsl5355555 to be sig- 
nificantly associated with the antibody response to NDV. 

Bioinformatic analysis of candidate regions 

As shown in Additional file 3: Table S2, there are 194 
genes or non-coding RNAs in the candidate regions of 
1-Mb windows (SNP position ± 0.5 Mb) surrounding 
each selected SNP (P<1.00xl0~ 4 ). These genes participate 
in at least 25 KEGG pathways, and comprise at least 108 
protein domains in the Interpro database of chickens. 
Most of the genes in the significantly (P<0.1) enriched 
pathways and protein domains are related to the immune 
system (Additional file 5: Table S3). Among these, CD2 
and CD80, which contain nearly significant CP<1.00xl0" 4 ) 
SNPs (rsl3904539, rsl3904683 and GGaluGA032608) for 
the antibody response to NDV, belong to the cell adhesion 
molecules pathway (gga04514), which involves the major 
histocompatibility complex (MHC) and T cell receptor 
genes. LOC418424, junction adhesion molecule 2, CD 200, 
CD80, ROBOl and ROB02 on GGA1, and Aggrecan on 
GGA10 all possess a conserved protein domain called the 
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immunoglobulin^ fold (IPR013783, Additional file 5: 
Table S3), suggesting that these proteins could be mem- 
bers of the complex gene network that modulates the 
immune response. The Cbl proto-oncogene, E3 ubiquitin 
protein ligase B, inter ieukin 10 receptor, beta, interferon 
gamma receptor 1 and interferon gamma receptor 2 genes, 
whose proteins participate in the Jak-STAT signaling path- 
way (gga04630), are also implicated in the antibody re- 
sponse to NDV. 

Discussion 

The host immune response to viruses is a complex 
process. The antibody response to NDV can be considered 
to be a quantitative trait under polygenic control, but 
with some QTLs [10,11]. In the current study, GWAS 



identified a region located approximately 100 Mb from the 
proximal end of GGA1, in which rsl5354805 was most 
significantly CP<1.26xlO~ 6 ) associated with, and accounted 
for 5% of the phenotypic variation of, the antibody response 
to NDV (Table 1). Moderately high heritability (0.41 ± 0.17) 
of the chicken antibody response to NDV was observed 
three decades ago [15]. Heritabilities of the primary anti- 
body response to NDV two weeks post-vaccination were 
0.27 ± 0.06 and 0.29 ± 0.05 in two Tanzania chicken eco- 
types {Kuchi and Tanzania Medium), respectively [16]. 
Lwelamira also reported recently that the heritability of the 
chicken primary antibody response to NDV at two weeks 
post-immunization was 0.22 ± 0.08 in the offspring of 
Kuchi chickens [17], a figure comparable to the 0.24 ± 0.08 
estimated here for the secondary antibody response. The 



Table 1 Genome-wide significant' SNPs for antibody response to the Newcastle disease virus 



SNP ID 


GGA 


Position (bp) 


Nearest gene 


SNP 


FA 


FAW 


FAR 


FAS 


P value 


R 2 


rs1 5354805 


1 


100399530 


ROB02 


CfT 


C 


0.756 


0.750 


1.000 


3.90x1 0" 7 


0.05 


rs1 5355555 


1 


100994585 


243 Kb upstream of ROB02 


A/G 


G 


0.758 


0.750 


1.000 


8.1 5x1 0" 7 


0.04 



Note: TA'=favorable allele with higher antibody response to the Newcastle disease virus, 'FAW'= favorable allele frequency of the F 2 population, TAR'= favorable 
allele frequency of the fast-growing Chinese yellow broiler sire line, TAS'= favorable allele frequency of the Huiyang beard chicken. ^Significant at the 5% level by 
the Bonferroni correction. 
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Figure 2 (See legend on next page.) 
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(See figure on previous page.) 

Figure 2 Genome scan for antibody response to the Newcastle disease virus on chicken {Gallus gallus) chromosome 1 with rsl 5354805 
and rsl 5355555 SNPs. The black line indicates genome-wide 5% significance, with a P-value threshold of 1.26x10" 6 Adjusted P values on a log 10 
scale for SNPs examined for their association with antibody response to the Newcastle disease virus. Structures, orientations and locations of the 
LOC769250, R0B01 and R0B02 genes on the NCBI Reference Sequence Build 2.1 are shown, together with pairwise LD estimates and haplotype 
blocks of the SNPs located within a 2-Mb region with rs15354805 and rsl 5355555 SNPs (red box). 



QTL linked to rsl5354805 therefore accounted for >20% 
(0.05/0.24) of the additive genetic variance of this trait. It is 
reasonable to consider that this QTL would be important 
for improving the ability of the immune response to NDV 
in chickens by use of marker-assisted selection. In addition, 
the effect of combining rsl5354805 on GGA1 and 
rsl4008095 on GGA10 is almost same as the sum of the 
effects of all SNPs with P values less than LOOxlO" 4 for 
the antibody response to NDV (7% vs. 8% of the pheno- 
typic variation) in this study. On that basis, there are prob- 
ably only two divergent QTLs for this trait on GGA1 and 
GGA10, while other SNP effects might result from LD. 

It should be noted that the region found in the present 
study to be significantly (P<1.26xl0" 6 ) associated with 
the antibody response to NDV differs from those found 
in earlier reports [10,11]. Yonash et al. reported that 
QTLs for the antibody response to NDV were located 
on GGA2 and GGA18 [11], whereas Biscarini et al. 
found 13 QTLs associated with the antibody response to 
NDV on GGA3, GGA4, GGA5, GGA9, GGA13, GGA16, 
GG A22 and GG AZ [10]. There are many possible rea- 
sons for the differing results among these studies. The 
population specificity of QTLs for the antibody response 
to NDV is probably the most important reason [18]. The 
current study was based on a broiler resource population, 
while the other two studies were based on nine different 
laying lines [10], and a chicken resource population built 
from the intercross of broilers with high and low antibody 
response to Escherichia coli vaccine [11]. In addition, 
Yonash et al. focused on the primary antibody response to 
NDV [11], while the current study examined the secon- 
dary antibody response to NDV. In the primary antibody 
response, the predominant class of antibody produced is 
immunoglobulin M, and in the secondary antibody re- 
sponse it is immunoglobulin Y; therefore, the ranking of 
birds on antibody response to the primary immunization 
may differ from that after the secondary boosting resulting 
in different QTLs being detected. In fact, the QTLs for 
the antibody response to NDV in this study may reflect 
the ability of the memory cell pool to respond to NDV. 
Interestingly, we observed that certain immunity- related 
genes, such as CD2, CD80 and CD200, on the detected 
QTLs (Additional file 3: Table S2). However, some 
detected QTLs in the current study were similar to those 
in two other studies, if the significance threshold set as 
P<0.05 (Figure 1). For example, rsl4520447 (GGA5: 
21181095 bp, P=0.0031), rsl4S20S28 (GGA5: 21329838 



bp, P=0.0194) and GGAluGA278144 (GGA5: 21355363 
bp, P=0.0194) were to the strongest of the QTLs detected 
by Biscarini et al. [10], and rsl0727400 (GGA18: 
4572183 bp, P=0.0336), rslS821136 (GGA18: 4587942 
bp, P=0.0171) and GGAluGA120237 (GGA18: 4594180 
bp, P=0.0200) were similar to the QTLs detected on 
GG Al 8 in the report of Yonash et al. [11]. 

The MHC plays an important role in animal immunity 
[19-22]. However, the present study found that the most 
reliable QTL linked to antibody response to NDV could 
be in the region that included rsl 5354805 and rsl 53 55 55 5 
on GGA1, whereas it did not detect any significant SNPs 
on GGA16, which harbors the chicken MHC. This may 
be because very few SNP markers on GGA16 could be 
used for QTL mapping in the current study (Additional 
file 2: Table SI). Further studies are necessary in order to 
detect QTLs for the antibody response to NDV on GGA16. 

The 2 Mb region of 99-101 Mb region containing 
three coding genes on GGA1 had the most significant 
effects for the secondary antibody response to NDV 
(Figure 2). Within this region, LOC769250 is a predicted 
gene with no clear function, and ROB 01 and ROB 02 
belong to the ROBO family, which are members of the 
immunoglobulin superfamily [23]. The encoded prod- 
ucts of both genes are integral membrane protein recep- 
tors for the SLIT-family, and have similar functions in 
axon guidance and migration of neuronal precursor cells 
[24-27]. Many molecules, such as Semaphorin 7A, 
Neurokinin A and Plexin-B2, participate in nervous sys- 
tem and immune system processes [28-32]. Perhaps 
ROBOl and ROB02 also affect the immune system. In 
fact, SLIT/ROBO signaling can influence the expression 
of Racl and CDC42 [33]. Racl plays a crucial role in T 
cell development [34] while CDC42 is essential for mem- 
ory T cell growth and differentiation [35,36]. Therefore, 
SLIT/ROBO signaling could contribute to adaptive im- 
munity through Racl or/and CDC42. Consistently, the 
focus of QTL mapping in this study was the difference of 
secondary antibody response to NDV that would be 
mainly related to activation, growth and differentiation of 
the memory T cells after the second immunization. In 
addition, SLIT/ROBO signaling plays important roles in 
leukocyte chemotaxis [37-39]. Therefore, ROBOl and 
ROB02 could well be candidate genes for modulating the 
antibody response to NDV. However, determining the 
functional role of these QTLs in the antibody response to 
NDV needs requires further studies. Studies of chicken 
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network pathways involving SLIT/ROBO signaling need 
to be replicated to validate the present findings. 

In addition to ROBOl and ROB02, some of the other 
genes that are close to the QTL region (P<L00xl0~ 4 , 
Additional file 3: Table S2) in the current study might 
modulate antibody responses to NDV, based on their 
biological function and networks. These genes could 
comprise a complex network (Additional file 6: Figure S3). 
On GG A10, interleukin 16, whose encoded protein is 
linked to ROB02 by three pathways (Additional file 6: 
Figure S3), participates in leukocyte chemo taxis and posi- 
tive regulation of responses to external stimuli [40], while 
B-cell CLL/lymphoma 2 - related protein Al is essential 
for lymphocyte activation and cell survival [41]. In fact, 
seven genes in the QTL region on GGA1 were enriched 
in the cell adhesion molecules pathway, which affects the 
immune system (Additional file 5: Table S3). Among these 
genes, CD2 and CD80 were initially considered to be can- 
didate genes. CD2 not only plays an important role in the 
generation of nanotubes between natural killer cells and 
their targets [42], but also provides signals to control the 
development and differentiation of T cells [43-45]. CD80 
is expressed on B cells, but also induces T cell proliferation 
and cytokine production [46-48]. In addition, CD200, 
which is about 180 Kb from rsl48S6616 (P=6.08xl0~ 6 ), is 
expressed on the surface of immune cells and affects the 
adaptive immune system by stimulating proliferation of T 
cells [49]. Interestingly, evidence shows that dendritic cells 
and T cells show a stronger response in the absence of 
CD200 inhibitory signaling [50]. In addition, there is a re- 
lationship between CD200 and influenza virus infection 
[51]. SAM domain, SH3 domain and nuclear localization 
signals 1 is another immune- related protein that partici- 
pates in B cell activation and differentiation [52]. CD247, 
IFNGR2 and IFNAR1 belong to the natural killer cell 
mediated cytotoxicity pathway (Additional file 5: Table S3). 
Obviously, a large number of genes related to innate and 
adaptive immunity are in the 89.7-111.0 Mb region of 
GGA1. This QTL region also overlaps the QTL for the 
antibody response to sheep red blood cells [53]. Thus 
this genomic region is of great importance for the 
chicken immune response and, possibly, disease resistance, 
and has been impacted by intense artificial selection (e.g., 
ZH p =-6.2 in ROB02) during chicken domestication [54]. 

Conclusions 

This study has exposed a novel QTL region associated 
with the antibody response to NDV. The QTL region 
was located approximately 100 Mb from the proximal 
end of GG Al. This region may play an important role in 
the chicken immune response. ROBOl and ROB02 were 
found to be important prospective candidate genes for 
modulating the antibody response to NDV in chickens; 
however, the roles of ROBOl and ROB02 require further 



study, such as exploring the effects of ROBOl and/or 
ROB02 silencing and overexpression in vitro and in vivo 
on immune gene networks. 

Methods 

Ethics statement 

The Animal Care Committee of the Institute of Animal 
Science, Guangdong Academy of Agricultural Sciences 
(Guangzhou, P. R. China) approved this study (approval 
number GAAS-IAS-2009-73). Animals involved in the 
current study were humanely euthanized, as necessary, 
to reduce any suffering. 

Experimental population 

The F 2 resource population was built from the full-sib 
intercross of two divergent lines (23 P and 51 F^. The first 
line was a fast-growing Chinese yellow broiler raised at 
Guangdong Wiz Agricultural Science & Technology Co. 
Ltd., which had undergone more than ten generations of 
selection for high growth rate and meat quality tailored to 
Chinese tastes. The second line was the Huiyang beard 
chicken, a Chinese local breed with a low growth rate and 
high meat quality tailored to Chinese tastes. The resis- 
tance of the second line is better than that of the first line 
(Data not shown). The F 2 population comprised 511 indi- 
viduals from six hatches. Before 40 days of age, each hatch 
was maintained in a group cage. From 40 days to 91 days, 
each chicken was reared in an individual cage with an 
individual trough. The house was equipped with water 
curtain systems and supplied with 24-hour light per day. 
A starter diet (200 g total protein and 2,900 kcal ME/kg) 
from hatch to 35 days of age and a grower diet (180 g total 
protein and 2,950 kcal ME/kg) from 43 days to 91 days of 
age were provided ad libitum. 

Phenotypic measurements 

The chickens were immunized with a commercial NDV 
live vaccine of the LaSota strain (Intervet International 
B.V., Boxmeer, Netherlands), using the standard dose 
given in the instructions of the vaccine, by eye drop at 
25 days of age and at 50 days of age, respectively. At 41 
days after the second immunization (91 days of age), 
serum samples were collected from the chickens. The 
antibody response to NDV was determined using an indir- 
ect ELISA and expressed as the S/P value of corresponding 
dilutions, according to the instruction of the commercial 
ELISA kit for testing the antibody response to NDV 
(BioCheck, Inc., Foster City, CA, USA). The ELISA kit 
measures NDV-specific immunoglobulin, Y but not im- 
munoglobulin A nor immunoglobulin M. 

Genotyping procedures 

Genomic DNA extraction from venous blood was per- 
formed using a phenol/chloroform method. For all 511 F 2 
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individuals, the quality and concentration of genomic 
DNA fulfilled the requirements of the Illumina Infinium 
SNP genotyping platform. Genotyping using the Illumina 
60K Chicken SNP Beadchip [12] was carried out at the 
Illumina-certified service provider, DNA LandMarks Inc., 
Saint-Jean-sur-Richelieu, Quebec, Canada. Quality control 
was assessed in GenomeStudio v2008.1 [55]. Six samples 
were excluded because they had more than 5% missing 
SNP genotypes. The final SNP set used in this study 
consisted of 39,833 SNP markers, meeting the following 
selection criteria: call frequency >95%, heterozygosity 
cluster intensity and separation value >0.4, and minor 
allele frequency >0.1. Information about the SNPs on each 
chromosome is summarized in Additional file 2: Table SI. 
Gene locations were based on the Ensembl Genome 
Browser [56] and the National Center for Biotechnology 
Information (NCBI) [23]. 

Data analysis 

Variance components were estimated using the average 
information restricted maximum likelihood algorithm 
[57] implemented by the DMU package [58]. The model 
was as follows: 

y = Xb + Za + e 

where y is the vector of observations of antibody response 
to NDV; b is the vector of fixed effects including sex (two 
levels) and hatch (six levels); a is the vector of animal addi- 
tive genetic effects; e is the vector of random residuals; X 
and Z are corresponding incidence matrices. The variance 
matrices of random effects were var(a) = G = A&1 and 
var(e) = R = \o 2 e , and the distributions for the random ef- 
fects were assumed to be: 

a~A/X0,G)e~A/ r (0,R), 

where A is the additive relationship matrix between 
individuals in the pedigree of all birds; o\ is animal addi- 
tive genetic variance; I is an identity matrix; and o 2 e is the 
residual variance. Heritability of antibody response to 

NDV was estimated as h 2 = 4 = and R 2 = 

thus, the ratio of the QTL effect compared with the addi- 
tive genetic variance was calculated as Q 2 = = , 
where (T 2 , (Tq TL and R 2 are the phenotypic variance, the 
QTL effect variance, and coefficient of correction for the 
regression equation with one or more QTL as indepen- 
dent variable and the antibody response to NDV as a 
dependent variable, which indicated the ratio of the QTL 
effect compared with the phenotypic variance of antibody 
response to NDV, respectively. 

A two-step GWAS procedure was performed in this 
study. First, the phenotypic data in the F 2 population 



were adjusted for sex, hatch, sire and dam effects by the 
following model: 

Yijki =ft + SEXi + Hj + S k + D kl + e ijk i 

where is the phenotypic value for the antibody re- 
sponse to NDV, is the overall mean, SEXi is the effect of 
the i th sex, Hj is the effect of the j th hatch, S k is the effect of 
the sire, D w is the effect of the 1 th dam within k** 1 sire, 
and e^ is the residual effect. The adjusted phenotypic value 
was calculated as ft + e^. 

Second, the adjusted antibody response to NDV, 
which followed a normal distribution (Additional file 1: 
Figure SI), was used to carry out genome- wide associ- 
ation analysis in PLINK using linear regression analyses 
without any covariates other than SNPs [59]. The thresh- 
old P value of the 5% Bonferroni genome- wide signifi- 
cance was 1.26xl0" 6 (0.05/39833), and the threshold P 
value for the "suggestive linkage" significance, which 
allowed one false positive effect in a genome-wide test, 
was 2.51xl0~ 6 (1/39833). Using the "gap" package in 
R v2.12.0 (www.r-project.org), a Manhattan plot was 
drawn to describe the P values of all SNPs. This study 
calculated the fraction of the phenotypic variance 
explained by the significantly associated SNPs according 
to a previously described method [60]. Haploview soft- 
ware [61] was used to calculate LD as an r-square value 
and build haplotype blocks under the four gamete rule 
(the 4 th gamete must be observed at frequency >0.01) in 
the candidate region associated with antibody response 
to NDV in the F2 population. 

Bioinformatics analysis 

Names of genes in the candidate regions with 1-Mb 
windows (SNP position ± 0.5 Mb) surrounding each 
significant SNP (P<L00xl0~ 4 ) were obtained from the 
Ensembl Genome Browser [56] and the NCBI [23]. A 
pathway and protein domain analysis was performed using 
the Database for Annotation, Visualization and Integrated 
Discovery (DAVID) [62,63]. A network of gene correlation 
was constructed by integrating the pathways from the 
human database of the Kyoto Encyclopedia of Genes and 
Genomes (KEGG) [64] (www.genome.jp/kegg/), the Gene 
Map Annotator and Pathway Profiler (GenMAPP) [65] 
(genmapp.org/) and the BioCarta (www.biocarta.com/ 
genes/index.asp). 

Additional files 

f \ 

Additional file 1: Figure SI. Distribution of the adjusted antibody 
response to the Newcastle disease virus in the F2 population. The data 
fitted the normal distribution (Shapiro-Wilk W test P=0.5527). 

Additional file 2: Table SI. SNP markers used in the genome-wide 
association study. 
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Additional file 3: Table S2. SNPs associated with antibody response to 
Newcastle disease virus (P<1.00x1CT 4 ). 

Additional file 4: Figure S2. Pattern of linkage disequilibrium on 
chicken (Gallus gallus) chromosome 1. 

Additional file 5: Table S3. Functional annotation of candidate genes 
by KEGG and Interpro analysis in chickens (P<0.1). 

Additional file 6: Figure S3. Network of gene correlation for genes in 
candidate regions related to antibody response to Newcastle disease 
virus (P<1.00x1CT 4 ). Yellow ellipses indicate genes. Green lines, red lines 
and blue lines indicate pathways based on GenMAPP, KEGG and BioCarta, 
respectively. The numbers indicate the number of pathways related to 
the two genes. 
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